function obj = callobj(C,Zp,Z,alpha,gamma,eta)
% Compute objective
numview = size(C,3);

obj = sum(C.*Zp,'all') + sum(Zp.^2,'all') * gamma + ...
    sum(sum(repmat(Z,1,1,numview)-Zp,[1,2]) .* alpha.^2) * eta;

end

